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There has never been a more exciting time in the overlapping areas of nuclear 
physics, particle physics and relativistic astrophysics than today. Orbiting ob- 
servatories such as the Hubble Space Telescope, Rossi X-ray Timing Explorer 
(RXTE), Chandra X-ray satellite, and the X-ray Multi Mirror Mission (XMM) 
have extended our vision tremendously, allowing us to see vistas with an unprece- 
dented clarity and angular resolution that previously were only imagined, enabling 
astrophysicists for the first time ever to perform detailed studies of large sam- 
ples of galactic and extragalactic objects. On the Earth, radio telescopes (e.g., 
Arecibo, Green Bank, Parkes, VLA) and instruments using adaptive optics and 
other revolutionary techniques have exceeded previous expectations of what can 
be accomplished from the ground. The gravitational wave detectors LIGO, LISA 
VIRGO, and Geo-600 are opening up a window for the detection of gravitational 
waves emitted from compact stellar objects such as neutron stars and black holes. 
Together with new experimental forefront facilities like ISAC, ORLaND and RIA, 
these detectors provide direct, quantitative physical insight into nucleosynthesis, 
supernova dynamics, accreting compact objects, cosmic-ray acceleration, and pair- 
production in high energy sources which reinforce the urgent need for a strong and 
continuous feedback from nuclear and particle theory and theoretical astrophysics. 
In my lectures, I shall concentrate on three selected topics, which range from 
the behavior of superdense stellar matter, to general relativistic stellar models, to 
strange quark stars and possible signals of quark matter in neutron stars. 

1 Introduction 

A forefront area of research, both experimental and theoretical, concerns the ex- 
ploration of the subatomic structure of superdense matter and the determination 
of the equation of state - that is, the relation .between pressure P, temperature 
T and density e - associated with such matter.EJ Knowing its properties is of key 
importance for our understanding of the physics of the early universe, its evolution 
in time to the present day, compact stars, various astrophysical phenomena, and 
laboratory physics. The high-temperature domain of the phase diagram of super- 
dense matter is probed by relativistic heavy-ion colliders. Complementary to this, 
neutron stars contain cold superdense matter permanently in their centers (cf. Fig. 
Q), which make them superb astrophysical laboratories for probiwi.thc low-density 
high-density domain of the phase diagram of superdense matter .Em 

Neutron stars are dense, neutron-packed remnants of massive stars that blew 
apart in supernova explosions. They are typically about 10 kilometers across and 
spin rapidly, often making several hundred rotations per second. The discovery 
rate of new rotating neutron stars, spotted as pulsars by radio telescopes is 
rather high. To date about 1400 pulsars are known. Depending on star mass and 
rotational frequency, gravity compresses the matter in the core regions of pulsars 
up to more than ten times (~ 1.5 GeV/fm 3 ) the density of ordinary atomic nuclei, 
thus providing a high-pressure environment in which numerous subatomic particle 
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Figure 1. Competing structures and novel phases of subatomic matter predicted by theory to 
make their appearance in the cores (R i; 8 km) of neutron starsfl 



processes plausibly compete with each other and novel phases of matter may exist. 
The most spectacular ones stretch from the generation of new baryonic particles 
(e.g., S, A, 5, A) to quark (it, d, s) deconfinement to the formation of Boson conden- 
sates (n~ , K ~, H-matter), as illustrated in Fig. |l|. There are theoretical suggestions 
of even more exotic processes inside pulsars, such as the formation of absolutely 
stable quark matter, a configuration of matter even more stable than the most sta- 
ble atomic nucleus, 56 Fe! In the latter event, pulsars would be largely composed 
of pure quark matter, eventually enveloped in thin nuclear crusts (bottom-left por- 
tion in Fig. Q). No matter which physical processes are actually realized inside 
neutron stars, each one leads to fingerprints, some more pronounced than others 
though, in the observable stellar quantities. Paired with the unprecedented wealth 
of new observational pulsar data, it seems to be within reach for the first time 
ever to seriously explore the subatomic structure of matter in the high-density low- 
temperature portion of its phase diagram from observed pulsar data. To this aim 
Einstein's field equation of relativistic gravity, 

C" = R" V --g tiU R = 8TTT' iU (e,P(e)), (1) 
is to be solved in combination with the latest theories of the subatomic structure 
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of matter a The latter follow according to the scheme 

-ir-^ =0 p(e) ' (2) 

where £({<fi}) denotes a given stellar matter lagrangianE In general, £ is a compli- 
cated functional of the numerous hadron and quark fields, collectively written as 
{4>}, that acquire finite amplitudes up to the highest densities reached in the cores of 
neutron stars. According to what has been said just above, plausible candidates for 
4> are thejiargecLstates of tli^n/SAliaryon octet, p, n, S, A, HQ, the charged states 
of the A Ella, 7T- H and K~ m3tMM3 mesons, as well as u, d, s quarksJl§Mo The 
conditions of chemical equilibrium and electric charge neutrality of stellar matter 
require the presence of leptons too, in which case <f> = e~ , p~ . Theories of super- 
dense matter enter Einstein's field equation (|l|) via the energy-momentum tensor 
T^ v , which contains the equation of state of the stellar matter, P(e). Because of 
the rather uncertain behavior of the matter at supernuclear densities, the models 
derived for the equation of state differ considerably with from each other. This has 
its origin in various sources such as: (1) the many-body technique used for the de- 
termination of the equation of state, (2) the model adopted for the nucleon-nucleon 
force, (3) assumptions about the fundamental building blocks of neutron star mat- 
ter, (4) the inclusion of boson condensates, and (5) the considerations of a possible 
phase transition of confined hadronic matter into deconfined quark matter .0113 

In general Eqs. (jl|) and (||) were to be solved simultaneously since the particles 
move in curved spacetime whose geometry, determined by Einstein's field equations, 
is coupled to the total mass energy e of the matter. In the case of neutron stars, 
however, the long-range gravitational force can be cleanly separated from the short- 
range nuclear force so that Eqs. (|l|) and (||) decouple from each other. This simplifies 
the study of compact stellar objects considerably. 



2 Models for the Equation of State of Superdense Neutron Star 
Matter 

2.1 Nonrelativistic Models 

For non-relativistic models, the starting point are phenomenological nucleon- 
nucleon interactions, denoted Vij, which fit the nucleon-nucleon scattering data 
and the deuteron properties. In order to achieve the correct binding energy of nu- 
clear matter at the empirical saturation density po (po = 0.15 nucleons/fm 3 which 
corresponds to a mass density of e = 2.5 x 10 14 g/cm 3 ), two-nucleon potentials are 
supplemented with three-nucleon interactions Vijk- The hamiltonian is then of the 
form 

W = E(^) V ? + E^ + E ( 3 ) 

i ^ ' i<j i<j<k 

A proven many-body method frequently adopted to solve the associateiljpany-body 
Schroedinger equation 7^1^ >= E\^ > is the variational approach EjEJ where a 
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Figure 2. Graphical illustration of AV14 + UVII, UV14+UVII, UV14+TNI, G300, and Gf 00 J 



variational trial function \& v > is constructed from a symmetrized product of two- 
body correlation operators (Fij) acting on an unperturbed ground-state, i.e., 



i<j 



1$ > 



where |$ > denotes the antisymmetrized Fermi-gas wave function, 



(4) 



(5) 



The correlation operator contains variational parameters which are varied to mini- 
mize the energy per baryon for a given density pitj'EII 



E v {p) 



< | > 



> En 



(6) 



As indicated, E v constitutes an upper bound to the ground-state energy Eq. The 
energy density e(p) and pressure P(p) are obtained from Eq. (^) through the fol- 
lowing manipulations: 



e(p) = p (E v (p) + m) 



P(p) = p*^-E v (p), 



(7) 



which can be combined to an equation of state of the form P(e) applicable to the 
stellar structure calculations. 

A representative collection of non-relativistic nuclear equations of state is listed 
in Tables [l] and ^ where the following abbreviations are used: N = pure neutron 
matter; NP = n, p, and leptons in chemical equilibrium; and K denotes the incom- 
pressibility (in McV) of symmetric nuclear matter at saturation density. Several of 
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Figure 3. Graphical illustration of equations of state HV, HFV, FP(V 14 +TNI), G 30 o, and Gg™ 1 . 



Table 1. Non-relativistic models for the equation of state of neutron star matter. 



EOS Properties (see text) 


R 


jferences 


WFF(UVi4+TNI) NP, K=261 
WFF(UVi4+UVII) NP, K=202 
WFF(AVi4+UVII) NP, K=209 
FP(Vi4+TNI) N, A'=240 
APR(Ai 8 +UIX) NP 
TF96 N, if =234 


i'i 

22 
22 
23 
24 
25 


S0i 



these equations of state are shown in Figs. || and |3| where the pressure is plotted as 
a function of total energy density (in units of the density of normal nuclear mat- 
ter, e = 140 MeV/fm 3 ). The Thomas-Fermi equation of state TF26 of Table [l] is 
based on the new Thomas-Fermi approach of Myers and SwiateckiO The effectpzr; 
interaction w 12 of this new approach consists of the Seyler-Blanchard potential E3, 
gener.alized by the addition of one momentum dependent and one density dependent 
term,Eil 

Y(r 12 ) 



Vl2 



Po 



'7- 



Pl2 



,2p 



)}■ 



(8) 



The upper (lower) sign in (|^) corresponds to nucleons with equal (unequal) isospin. 
The quantities kp 0) To (= k 2 Fo /2m), and po are the Fermi momentum, the Fermi 
energy and the particle density of symmetric nuclear matter. The potential's radial 
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Table 2. Nuclear matter properties of the equations of state compiled in Table 



EOS 


E/A 


Po 


K 


M* 


Q-sy 




(MeV) 


(fm" 3 ) 


(MeV) 


(MeV) 


(MeV) 


WFF(UVi4+TNI) 


-16.6 


0.157 


261 


0.65 


30.8 


WFF(UVi4+UVII) 


-11.5 


0.175 


202 


0.79 


29.3 


WFF(AVi4+UVII) 


-12.4 


0.194 


209 


0.66 


27.6 


FP(Vi4+TNI) 


-16.00 


0.159 


240 


0.64 




APR(Ai 8 +UIX) 


-16.00 


0.16 








TF96 


-16.04 


0.161 


234 




32.0 



dependence is describe by the normalized Yukawa interaction 



i p—rii/a 

Y(r 12 ) = — . (9) 



Its strength depends both on the magnitude of the particles' relative momentum, 
P12, and on an average of the densities at the locations of the particles. The param- 
eters £ and £ were introduced in order to achieve better agreement with asymmetric 
nuclear systems, and the behavior of the optical potential is improved by the term 
a(2p/p ) 2 / 3 - Here the average density is defined by p 2 / 3 = (p^ 3 + p^ 3 )/2, where 
pi and pi are the relevant densities of the interacting particle (neutron or protons) 
at points 1 and 2. The potential's seven free parameters a, (3, 7, er, £, £, a are 
adjusted to the properties of finite nuclei, the parameters of the mass formula, and 
the behavior of the optical potentialO The nuclear matter properties at satura- 
tion obtained for TF96 are summarized in Table |^, where the listed quantities are: 
binding energy of normal nuclear matter at saturation density, E/A; compression 
modulus, K; effective nucleon mass, M* (= m* /m where m denotes the free nu- 
cleon mass); and the asymmetry energy, a sy . The new Thomas-Fermi force has the 
advantage over the standard Seyler-Blanchard interaction to not only reproduce the 
ground-state properties of finite nuclei and infinite symmetric nuclear matter but 
also the optical potential and, as revealed by stellar structure calculations E3, the 
properties of neutron stars as well. These features render the new Thomas-Fermi 
model very attractive for investigations of the properties of dense nuclear matter. 



2.2 Relativistic Fieldtheoretical Models 

Relativistic, field-theoretical theories of dense nuclear matter and finite nuclei have 
enjoyed a renaissance in recent years, and they have the virtue of describing nuclear 
matter at saturation, many features of finite miclei, both spherical and deformed, 
and they extrapolate causally to high density.EJ The starting point of relativistic 
theories of superdense neutron star matter is a lagrangian of the following typejj 



C(x) 



E 



C° B (x) 



(10) 



B=p,n,£±.°,A,E:o.-,A + + - + . .- 
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Table 3. Masses, quantum numbers (spin Jb, isospin Ib, strangeness Sb, hypercharge Yb, third 
component of isospin /3s), and electric charges (gs) of baryons. 



Baryon (B) 


m B (MeV) 


Jb 


Ib 


Sb 


Y B 


I3B 


QB 


n 


939.6 


1/2 


1/2 





1 


-1/2 





P 


938.3 


1/2 


1/2 





1 


1/2 


1 


£+ 


1189 


1/2 


1 


-1 





1 


1 


s° 


1193 


1/2 


1 


-1 











sr 


1197 


1/2 


1 


-1 





-1 


-1 


A 


1116 


1/2 





-1 











S° 


1315 


1/2 


1/2 


-2 


-1 


1/2 







1321 


1/2 


1/2 


-2 


-1 


-1/2 


-1 


A++ 


1232 


3/2 


3/2 





1 


3/2 


2 


A+ 


1232 


3/2 


3/2 





1 


1/2 


1 


A 


1232 


3/2 


3/2 





1 


-1/2 





A- 


1232 


3/2 


3/2 





1 


-3/2 


-1 



M=a,oj,n,g,r),8,4> [ B=p,n,...,A++.+> .- J L=p- .41- 

The subscript B runs over all baryon species that may become populated in dense 
neutron star matter (Table ||) . Thc_nuclear forces are mediated by a collection of 
scalar, vector, and isovector mesons Ej which is compiled in Table [|. The equations 
of motion that follow from Eq. ( [To|) for the baryon fields are given by 



(ij^dft - m B ) ipB (x) = 9oB <t(x) 4>b (x) 

A. 

IB 



fpB 

Am B 



+ {g P B r r-p^x) + a ^ t-G^(x)} i> B (x) 

+ hi ^ ( d T . n (x)) i/> B (x) , (11) 

where the standard field-theoretical notation of Ref. @ is used. The meson fields 
in (|ll]) obey the following field equationsJ3 

(d M d M + ml) <j(x) = - ^.gcrs 4>b{x)iI}b{x) - m N In 9oN {g a Nv{x)) 

B 

- CN 9<jN (gaNcr(x) f , (12) 
d^F^x) +m*<j l/ (x) = ^^B ^b{x)^^ b {x) 

B 

-^d^{i, B {x)a^ B (x))) 1 (13) 
(d^ + ml) tt(x) =J2—d» (4> B {x) 75 7m r ^ B {x)) , (14) 

B 
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Table 4. Mesons and their quantum numbers 
and mass mjif of meson M. 



s 



The entries are: spin Jm, parity tt, isospin Im, 



Meson (M) 


T7T 


Im 


Coupling 


Mass (MeV) 


cr 


0+ 





scalar 


550 


w 


1- 





vector 


783 


tt* 


Q- 


1 


pseudovector 


140 


^° 


o- 


1 


pseudovector 


135 


P 


1" 


1 


vector 


769 


?7 


o- 





pseudovector 


549 


<5 


0+ 


1 


scalar 


983 




1- 





vector 


1020 




o- 


1/2 


pseudovector 


494 


K- 


o- 


1/2 


pseudovector 


494 



d^G^x) + m 2 p p u (x) = ^2^g P B 4>b(x)tj u 4>b(x) 

B 



fpB 



^(feWrvfefa:))}, (15) 



2m B 

with the field tensors F^ v and denned as F» v = d^ui" - d v uj» and = 
Qfi pv _ qv Equations Jll| ) through ([lq ) are to be solved simultaneously in 
combination with the neutron star matter constraints of electric charge neutrality 
and chemical equilibrium. The former can be written as 



J2<1B (2J B + 1) 



Pf,b 
6tt 2 



E - QM&(» M -m M ) = 0, (16) 

A— e./i 



while chemical {(3) equilibrium is expressed as 

li B = /i™ - q B n e , 



(17) 



where /i™ and /i e denote the chemical potentials of neutrons and electrons whose 
knowledge is sufficient to determine the chemical potential of any other baryon 
present in the system. The only rnesons that may pl a u s i b ly condense (density pm) 
in neutron star matter are the n~ Bo or the K~ E9l3E30. Equations (|16|) and ([17]) 
leads to constraints on the Fermi momenta of baryons and leptons, kp B and kp L 
respectively. The leptons in neutron star matter can be treated as free particles, 
which obey the following Dirac equation, 



(ry^ -m L ) ip L (x)=0, 



(18) 



with L = e~ , ji~ . An elegant mathematical framework that enables one to solve the 
equations of motion ( pT[ ) through (|l^) is the relativistic Greens function method, 
which is outlined next. The starting point is the Martin-Schwinger hierarchy of 
coupled Greens functions.ES In the lowest order, the Martin-Schwinger hierarchy 
can be truncated by factorizing the four-point Greens function 172(1, 2; 1'2'), where 
primed (unprimcd) arguments refer to ingoing (outgoing) fcrmions, into a product 



paper: submitted to World Scientific on February 5, 2008 




[ 2/ 






X 




X 


2' 


V 





g 2 (1,2;1\2') g^l.T) g 1 (2,2') 



2' 



Figure 4. Factorization of four-point baryon Green function 32 into antisymmetrized products of 
two-point baryon Green functions, g\ xgi. Direct (Hartree) and exchange (Fock) contribution are 
shown. This factorization scheme truncates the many-body equations at the Hartree-Fock level. 



of two-point Greens functions 1'), as schematically shown in Fig. |j. This leads 
to the well-known relativistic Hartree (i.e., mean-field) and Hartree-Fock approx- 
imations. The T-matrix approximation, also known as ladder (A) approximation, 
goes beyond Hartree and Hartree-Fock by truncating the Martin-Schwinger hier- 
archy by factorizing the six-point Greens function ^3(123; 1'2'3') into products of 
four-point and two-point functions by means of which dynamical two-particle cor- 
relations in matter, which are connected to a two-body boson-exchange potential v, 
are taken into account. The T matrix (effective two-particle interaction in matter) 
is computed from an integral equation of the following type, u 

T = v -v cx + J v AT , (19) 

with A the baryon-baryon propagator in intermediate scattering states. In the 
framework of the so-called A 00 approximation, the baryons in intermediate states 
propagate as free particles, leading to a A w A 00 = iglg®. This is in sharp con- 
trast to the relativistic Brucckner-Hartree-Fock (RBHF) approximation where the 
intermediate baryons are coupled to the nuclear background and, therefor e are con- 
strained to scattering states outside their respective Fermi seas.oEJo The basic 
input quantity in Eq. ( [l9| ) is the nucleon-nucleon interaction in free space. Rep- 
resentative models are, for instance, the Bonn meson-exchange model and the 
boson- exchange potentials of Brockmann and MachleidtHa A characteristic fea- 
ture of these potentials is that their parameters are adjusted to the relativistic 
two-nucleon scattering data and the properties of the deuteron. In this respect a 
parameter-free treatment of the many-body problem is achieved. The Born approx- 
imation to the T-matrix, which merely sums numerous boson-exchange potentials, 

< 12 I v I 1'2' > = Y, r "' r 22< 4* . (20) 

neglects dynamical correlations among baryons. It is this approximation which 
reduces the ladder approximation to the relativistic Hartree and Hartree-Fock 
approximations^ The symbols T M in Eq. (^0j) stands for the various meson-nucleon 
vertices, and A M denotes the free meson propagator of a meson of type M whose 
explicit forms are given Ref. The baryon self-energy, or mass operator, is obtained 
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Figure 5. Graphical representation of g^ B T, B g B of Dyson's equation (|22|). The Fermi sea of each 
baryon is filled up to the highest single-particle energy state, . The infinite Dirac sea is filled 
with antibaryons, which too modify the motion of the baryons. The Fermi an d Rim e-sea graphs 
are referred to as medium and vacuum polarization contributions, respectively. P !! 43 ! 



tr [T BB ' gf) -T BB ' gf 



(21) 



from the T matrix according tola 

B'=p,n,S±- ,A,S .-,A+ + . + .°.' 

where B' sums all the charged baryon states whose thresholds are reached in neutron 
star matter treated within this framework. The baryon Greens functions g B in Eq. 
(Uw are given as solutions of Dyson's equation, 

9? - 9l B + 9 r^ B {{g B '}) 9 B , (22) 

which terminates the set of equations to be solved self-consistently. The diagram- 
matic representation of the second term of Dyson's equation is shown in Fig. ||. It 
is this term which corrects the free Greens function for medium effects originating 
from the presence of Fermi seas of filled baryons. The equation of state follows 
from the stress-energy density tensor, T^, of the system asS 

E(p) = < Too > IP ~ m , where (23) 

T, v {x) = £ ^ - 9, v L{x). (24) 

X—B,L 

The pressure is obtained from E(p) as in Eq. (0). 

A broad collection of relativistic models for the equation of state of superdense 
neutron star matter is compiled in Table ||. The specific properties of these equa- 
tions of state are described in Table ||, where the following abbreviations are used: N 
= pure neutron; NP = n, p, leptons; it = pion condensation; H = composed of n, p, 
hyperons (S^ , A, S 0,_ ), and leptons; A = Ai 2 32-resonance; K~ = condensate of 
negatively charged kaons; and K = incompressibility (in MeV). All equations of 
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Table 5. Relativistic, field-theoretical models for the equation of state of neutron star matter. 



EOS 



pUCMi 

"225 

pDCM2 

"265 

pK240 

"M78 

G30O 

r<K~ 

^300 



Properties (see text) References 



if 



G7T 
200 

"300 

HV 

^Bonn 

HFV 

A 00 . 
il HEA 
A BHF . 
^BroB 



-HV 

HFV 
HFV 



H, K 
H, K 
H, K 
H, K 
R,K~ 

H, 7T, 
H, 7T, 

H, X: 

H, K 
H, A, 
H, A, 
H,A, 



=225 

=265 

=240 

=300 

, X=300 

if=200 

£T=300 

=285 

=186 

X=376 

#=115 

#=249 



r 

IS 



Table 6. Nuclear matter properties of the equations of state compiled in Table 



EOS 




E/A 


Po 


K 


M* 


G-sy 






(MeV) 


(fm- 3 ) 


(MeV) 


(MeV) 


(MeV) 


pUCMl 
"225 




-16.0 


0.16 


225 


0.796 


32.5 


nDCM2 
"265 




-16.0 


0.16 


265 


0.796 


32.5 


pK240 
"M78 




-16.3 


0.153 


240 


0.78 


32.5 


G30O 




-16.3 


0.153 


300 


0.78 


32.5 


r\K~ 
"300 




-16.3 


0.153 


300 


0.78 


32.5 


f-17T 

"200 




-15.95 


0.145 


200 


0.8 


36.8 


"300 




-16.3 


0.153 


300 


0.78 


32.5 


HV 




-15.98 


0.145 


285 


0.77 


36.8 


A-Bonn "+ 


HV 


-11.9 


0.134 


186 


0.79 


34 


HFV 




-15.54 


0.159 


376 


0.62 


30 


^HEA + 


HFV 


-8.7 


0.132 


115 


0.82 


29 


ABHF , 
^BroB ' 


-HFV 


-15.73 


0.172 


249 


0.73 


34.3 



state of this collection account for neutron star matter in full baryonic equilibrium. 
A few selected equations of state of this collection are shown in Fig. |[ An inherent 
feature of relativistic equations of state is that they do not violate causality, i.e., 
the velocity of sound given by v s = \J dP jde is smaller than the velocity of light at 
all densities. This is not the case for most non-relativistic models for the equation 
of state. Among the latter, WFF(UV 14 + TNI), APR(A 18 + UIX), and the TF96 
models for the equation of state do not violate causality up to densities relevant for 
the construction of models of neutron stars from them. 
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2.3 Baryon-Lepton Composition of Neutron Star Matter 

At densities lower than the density of normal nuclear matter, eo (= 140 MeV/fm 3 ), 
neutron star matter consists of only p, n, e _ , and /i~ whose densities are determined 



by the matter equations (16) through (22). At densities larger than eo the more 



massive baryon states E^ , A, S 0, ~, A + +^'° become populated since their chemical 
potentials, [i B , become larger than the lowest-lying energy eigenstates 0J B (p = 0) 
of these particles, 

ix B =^ n -q Bf i e > uj B (p = 0). (25) 
In relativistic field theory, the energy eigenstates are given byi 

co B (p) = J$(w B (p),p)+I 3B X b 3 (lo b (p), P ) + [m B + Ef (cu B (p),p)] , (26) 
while in the non-relativistic case one has 

Lo B (p) = -^-p 2 + Z B (p), (27) 
2m B 

where E B denotes the non-relativistic one-particle potential felt by a baryon in 
matter. Since the chemical potentials pL n and /i e are positive and the self-energies 
Yi B 3 are negative (E^ is proportional to the difference of proton and neutron den- 
sities, which is negative for neutron star matter because of its neutron excess), it 
follows from Eq. ( p^ ) that (i) negatively charged baryons are charge-favored, and 
(ii) baryons having_the opposite isospin orientation as the neutron (i.e., I^b > 0) 
are isospin-favoredo The physical reason behind item (i) is that negatively charged 
baryons replace electrons with high Fermi momenta, while (ii) is making the matter 
more isopin symmetric. Both items enable neutron star matter to settle down in an 
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parameter set (see Ref. □). The hyperons are universally coupled hypcrons and A's.| 
non-universally coupled, the A's universally. 



energetically more favorable state. Figures |6| through |9| illustrate possible particle 
compositions of neutron star matter computed for different many-body techniques 
as well as different nuclear forces. One sees that in relativistic Hartree (HV) the 
A hyperon, which remains unaffected by the effects (i) and (ii) since its electric 
charge q\ = and the third component of isospin /3A = 0, has the lowest threshold 
of all hyperons. In the case of relativistic Hartree-Fock (HFV), the charge favored 
but isospin unfavored S~ possesses the lowest threshold, and the A makes its ap- 
pearance too. The latter is also the case for relativistic Brueckner-Hartree-Fock 
type calculations u The important net effect of all these new degrees of freedom is 
a softening of the equation of state at supernuclear densities, which reduces the 
maximum mass of a neutron star. 



2.4 H-dibaryons 

A novel particle that may make its appearance in the center of a neutron star 
is the H-dibaryon, a doubly strange six-quark composite with spin and isospin 
zero, and baryon number twoJla Since its first prediction in f977, the H-dibaryon 
has been the subject of many theoretical and experimental studies as a possible 
candidate for a strongly bound exotic state. In neutron stars, which may contain 
a significant fraction of A hyperons (Figs. ^| through ||), the A's could combine to 
form H-dibaryons, which could give way to the formation of H-matter at densities 
somewhere between 3 £q and 6 eo depending on the in-medium properties of the 
H-dibaryon.E2lE2Jc3 This could cause a significantpSoftening of the equation of state 
as shown in Fig. |l0| (curves labeled Uh)- In Ref. EJ it was shown that H-dibaryons 
with a vacuum mass of about 2.2 GeV and a moderately attractive potential in 
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Figure 10. Models for the equation of state of neutron star matter la 

the medium of about Uh = —30 MeV could go into a Bose condensate in the 
cores of neutron stars if the limiting star mass is about that of the Hulse- Taylor 
pulsar PSR 1913+16, M = 1.444 M@. Conversely, if the medium potential were 
moderately repulsive, around Uh = +30 MeV, the formation of H-dibaryons may 
only take place in heavier neutron stars of mass M ^ 1.6 Mq. If indeed formed, 
however, H-matter may not remain dormant in neutron stars but, because of its 
instability against compresaipii-cpuld trigger the conversion of neutron stars into 
hypothetical strange stars .ljSEJ 

2.5 Condensation of K~ Mesons 
Once the reaction 

e-^K~+v (28) 

becomes possible in a neutron star, it becomes energetically advantageous for the 
star to replace the fermionic electrons with the bosonic K~ mesons. Whether or 
not this actually happens depends on how quickly thp_ftT~ drops with density in 
dense matter (Fig. |Il|) . An analysis of the KaoS data shows that the attraction 
from nuclear matter may bring the K~ mass down to m* K _ ~ 200 MeV at p ~ 3 pq. 
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Figure 11. Effective kaon mass in nuclear and neutron star matter 
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For neutron-rich matter, the K mass may drop asl 



m* K -(p)~m K - (1-0.2-^) , 

Po 



(29) 



with rriK = 495 MeV the K~ vacuum mass. Combining this relation with the 
outcome for shown in Fig. [ll] shows that the threshold condition for the onset 
of K~ condensation, fx e = m* K would be fulfilled at densities p ^ 3po> which are 
easily reached in neutron stars. 

2.6 Quark Deconfinement 

The phase transition between confined hadronic matter and quark matter is char- 
acterized by the conservation of baryon charge and electric charge. The Gibbs 
condition for phase equilibrium then is that the two associated chemical poten- 
tials /i" and p e , corresponding to baryon and electric charge conservation, and the 
pressure in the two phases be equal, 



(30) 



Here Pr denoted the pressure of hadronic matter computed for a hadronic La- 
grangian £({(/>}) as givenJn Eq. ©. The pressure of quark matter, Pq, is obtain- 
able from the bag model.liZI The quark chemical potentials f/, u ,/j, d , fi s are related to 
the baryon and charge chemical potentials as 



3^ 



3^ 



3" 



3 M 



(31) 



Equation (30) is to be supplemented with the two global relations for conservation 
of baryon charge and electric charge within an unknown volume V containing A 
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Figure 12. Sample composition of chemically equilibrated neutron (quark-hybrid) star matter .| 




baryons.EJ The first one is given by 

P = y = (1 - X) PhOA /A T) + xPqW, /A T) , (32) 

where x = ^q/^ denotes the volume proportion of quark matter, Vq, in the un- 
known volume V, and pi (i = H, Q, L) is the particle number density. Global 
neutrality of electric charge within the volume V can be written as 

= | = (1 - x) ?h(m", /A T) + X «q(/i n , /A T) + q L , (33) 

with <7i the electric charge densities of hadrons, quarks, and leptons. For a given 
temperature T, Eqs. (^) through ( |33|) serve to determine the two independent 
chemical potentials and the volume V for a specified volume fraction x of the quark 
phase in equilibrium with the hadronic phase. After completion Vq is obtained as 
Vq = xV- Through Eqs. ( |30"| ) to (|3^), the chemical potentials obviously depend on 
the proportion x of the phases in equilibrium, and hence so also all properties that 
depend on them, the energy densities, baryon and charge densities of each phase, 



Table 7. Models for the equation of state of neutron star matter accounting for quark deconfine- 
ment. 



EOS Properties (see text) 


11 




rence 


G^fso 11 H,Q, if =225, i? 1/4 = 180 
Gg™ 2 H,Q, if =265, B 1 / 4 = 180 
G| 2 8 4 ° H,Q, if =240, BV4 = 180 
Gbiso h >Q> #=300, B 1 ' 4 = 180 


48 
43 


63 
63 
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Table 8. Nuclear matter properties of the equations of state compiled in Table 



EOS 


E/A 


Po 


K 


M* 


d S y 




(MeV) 


(£m- 3 ) 


(MeV) 


(MeV) 


(MeV) 


nDCMI 
"B180 


-16.0 


0.16 


225 


0.796 


32.5 


piDCM2 
"B180 


-16.0 


0.16 


265 


0.796 


32.5 


pK240 


-16.3 


0.153 


240 


0.78 


32.5 


pK300 
'-'BISO 


-16.3 


0.153 


300 


0.70 


32.5 



and the common pressure. For the mixed phase, the volume proportion of quark 
matter varies from < x < 1? an d the energy density is the linear combination of 
the two phases 

e = (1 - X) e H (M", M e , M, T) + X e Q ( M " , //, T) . (34) 

By solving the models of confined and dcconfincd phases in both pure phases and 
in the mixed phase, we can compute the baryon, lepton and quark populations 
in neutron star matter from Eqs. ( |30| ) through (|33"|). One such sample outcome 
is shown in Fig. [l^. Three features emerge immediately from this population. 
Firstly, one sees that the transition from pure hadronic matter to the mixed phase 
occurs at rather low density of about 3 po- Depending on the bag constant and the 
underlying nuclear many-body— approximation, threshold values even as small as 
about 2 p can be obtainedHS Secondly, the lepton density saturates as soon as 
deconfmed quark matter is generated, since charge neutrality can be achieved more 
economically among the baryon-charge carrying particles themselves. Thirdly, the 
presence of quark matter enables the hadronic regions of the mixed phase to arrange 
themselves to be more isospin symmetric than in the pure phase by transferring 
charge to the quark phase in equilibrium with it. Symmetry energy will be lowered 
thereby at only a small cost in rearranging the quark Fermi surfaces. Thus the 
mixed phase region consists of positively charged nuclear matter and negatively 
charged quark matter. Models for the equation of state of neutron star matter 
which account for quark deconfinement are listed in Tables ^ and |[ A graphical 
illustration of these equations of state is shown in Fig. [H3. 



3 Observed Neutron Star Properties 

Orbiting observatories such as the Hubble Space Telescope, the Chandra X-ray 
satellite, and the X-ray Multi Mirror Mission have extended our vision tremen- 
dously, allowing us to look at neutron stars with an unprecedented clarity and 
angular resolution that previously were only imagined, The properties of such ob- 
jects such as masses, rotational frequencies, radii, moments of inertia, redshifts, and 
temperatures are known to be sensitive to the adopted microscopic wpdel for the 
nucleon-nucleon interaction and thus to the nuclear equation of state.ErE-3 Hence by 
means of comparing theoretically determined values for these quantities with ob- 
served ones one may be able to constrain the behavior of superdense matter from 
the (unprecedented wealth of new) observational pulsar data. In the following we 
briefly summarize several important star properties. 
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Figure 13. Models for the equations of state of neutron star matter. The shai 
deconfinement taken into account in models G*l°°, Gg™ 2 , G^40 qDCMI |6 
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3.1 Masses 

The gravitational mass is of special importance since it can be inferred directly 
from observations of X-ray binaries and binary radio pulsar systems. The masses 
of the latter appear to be centered about 1.35 M Q .L3 The mass of the Hulse- Taylor 
radio pulsar PSR 1913+16 H, given by 1.444 ± 0.003) Mq, is somewhat higher 
than this value. Examples of neutron stars in X-ray binaries are Vela X-1 and 

!3@ and 



7+0.23 
-0.17 



Mr. 



the burster Cygnus X-2. Their respective masses are M = 1.87_ 
M = (1.8 ±0.4) Mq 13. Indications for the possible existence of very heavy neutron 
stars, with masses around 2 Mq (e.g., neutron star 4U 1636-536), may also come 
from the observation of quasi-periodic oscillations in luminosity in low-mass X- 
ray binaries.l-9 These mass values, which constitute upper boundaries on the star 
mass, are larger than the typical 1.35 Mq masses found in neutron star binaries, 
presumably due to accreted matter. 



3.2 Rotational Frequencies of Fast Pulsars 

The rotational periods of fast pulsars provide constraints on the equation of state 
when combined with the mass constraint. The fastest pulsars observed so far have 
rotational periods of 1.6 ms, which corresponds to a rotational frequency of 620 
Hz. The successful model for the nuclear equation of state must therefore account 
for rotational neutron star periods of at least P = 1.6 ms as well as for the masses 
discussed in Sect. 



3.1 



paper: submitted to World Scientific on February 5, 2008 



18 



Observationally, the searches for rapidly rotating radio pulsars are biased, be- 
ing least sensitive to short periods. Therefore the world's data of radip pulsars are 
unlikely to represent the true underlying population of fast pulsars £il Most of the 
large surveys have had very poor sensitivity to millisecond pulsars with rotational 
periods below about 4 ms, thus presumably distorting the statistics on pulsar pe- 
riods. The present cutoff in short periods at about 1 ms is therefore possibly only 
an artifact of the search sensitivity. 



3.3 Radii 

Direct radius determinations for neutron stars do not exist. However, combinations 
of data of 10 well-observed X-ray bursters with special theoretical assumptions lead 
Van Paradijs O to the conclusion that the emitting surface has a radius of about 8.5 
km. This value, however, may be underestimated by a factor of two.E3 Fujimoto and 
Taam derived from the observational data of the X-ray burst source MXB 1636- 
536, under somewhat uncertain-theoretical assumptions, a neutron star mass and 
radius of 1.45 M© and 10.3 km.u An error analysis lead them to predicting mass 
and radius ranges of 1.28 to 1.65 M Q and 9.1 to 11.3 Ism, respectively Very recently 
the nearby, isolated neutron star RXJ 185635-3754 LaO has attracted particular 
attention because of its interpretation as a possible strange quark star.lZ3 (The 
physics of strange quark stars will be discussed in Sect. ^.) The observed X-ray flux 
and temperature from this object, if originating from a uniform blackbody, appears 
to correspond to an effective stellar radius significantly smaller than ~ 10 km, the 
canonical neutron star radius. This follows from the relation 

/ o M \ _1 / 2 n 

R 00 = R[1 « 7 km , (35) 

\ R J 120 pc y ' 

where R is the neutron star radius, M its mass, and D the distance to the star. One 
major uncertainty in (^5|) is buried in the numerical coefficient, which could increase 
by 50% or more if the star's surface temperature is not uniform. Moreover, the 
existence of an atmosphere on this neutron star could alter the coefficient probably 
even more strongly.LS The other major uncertainty in obtaining a radius is -the star's 
distance which, however, has now been determined to be D = 117 ± 12 pc LJ, nearly 
double the originally published distance. Taking these considerations into accpimt, 
the interpretation of RXJ 185635-3754 as a quark star appears questionable. E2l 



3.4 Moment of Inertia 

Another global neutron star property is the moment of inertia, /. Early estimates 
of the energy- loss rate from pulsars spanned a wide range of J, i.e., 7 x 10 43 < I < 
7 x 10 44 g/cm 3 £3 From the luminosity of the Crab nebula (~ 2 — 4 x 10 38 erg/sec), 
several authors have found a low£r_bpund on the moment of inertia of the pulsar 
given by I> 4 - 8 x 10 44 g cm 2 BOB 
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3.5 Redshift 

Next we mention the neutron star redshift, z, which is given by 

/ 2M\-i/2 , 
-(l-^) -1, (36) 

with M and R the star's mass and radius. Liang has considered the neutron star 
redshift data base provided by measurements of 7-ray burst redshiftcd annihilation 
lines in the range 300 — 511 keV.EJ These bursts have widely been interpreted as 
gravitationally redshifted 511 keV pair annihilation lines from the surfaces of 
neutron stars. From this he showed that there is tentative evidence to support a 
neutron star redshift range of 0.2 < z < 0.5, with the highest concentration in 
the narrower range 0.25 < z < 0.35. A particular role plays the source of the 
1979 March 5 7-ray burst source, which has been identified with SNR N49 by its 
position. From the interpretation of its emission, which has a peak at ~ 430 keV, 
as the 511 keV annihilation line the resulting gravitational redshift has a value 
of z = 0.23 ±0.05. 

3. 6 Magnetic fields 

Neutron stars are highly magnetized objects.!! There is a general plausibility ar- 
gument to the effect that, if internal magnetic flux in a star is conserved owing 
to high conductivity, then the surface magnetic field increases as R~ 2 where R is 
the stellar radius. Thus, if the sun were to collapse to neutron star dimensions 
(roughly a factor of 7 x 10 4 in radius), one would magnify the general 1 G so- 
lar field to about 5 x 10 9 G. This value, while being in good agreement with the 
magnetic fielcLstrength of _B ~ 4 — 6 x 10 11 G estimated for the X-ray pulsar 
1 E 2259+586 E3, the central source in the supernova remnant G 109.1-1.0, is a bit 
shy of 10 12 G inferred from line features in the pulsed hard X-ray spectrum of 
Her X-l and 4U 0115-63.e3e3H Which physical mechanisms might be important 
to create the extraordinarily strong magnetic fields of neutron stars is an open is- 
sue. One possible mechanism may be differential rotation of the stellar core, as 
proposed for planetary magnetic fields. Besides that, ferromagnetism has been ex- 
amined as a possible magnetic field source, while differential rotation Jbetween say 
superfluid protons and normal electrons has been another suggestion. EJ Effects of 
strong magnetic fields on the structure of compact stars were studied, for instance, 
in Refs. ljO. Particularly striking are so-called magnetars, which are interpreted as 
neutron stars with superstrong surface magnetic fields on the order of ~ 10 15 G.E3 
This translates to interior magnetic fields of up to ~ 10 18 G, which is about the 
highest possible value a star can sustain. 

3. 7 Cooling Data 

The detection of thermal photons from a stellar surface via X-ray observatories 
like ASCA and ROSAT serves as the principal window on the properties of a star. 
The surface temperatures of stars are derivable from the measured flux and spec- 
trum. The cooling rate of a hot, young neutron star is primarily dependent at 
early times (first several thousand years) on the neutrino emissivity of the core's 
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Table 9. Overview of neutrino emitting processes relevant for neutron star cooling. I 



Name 


Particle 


processes 




Efficiency 


Modified Urea 


n + 


- n — 


> 71 + 


-P + 


e 


4 t^e 


slow 




n + 


■P + 


e~ - 


-> nH 


- 71 


4 V e 




Direct Urea 


n - 




- e~ 


4 P e 






fast 




P + 


e~ - 


-> n, 


4 1/ e 








Quark modified Urea 


d + 


a! 


e~ - 




-d- 


+ f e 


slow 






n + 


e _ - 


— > it - 


hd 


+ 






d + 


11 + 


e~ - 




- s - 








11 + 


n + 


eT - 


-*u- 


hs- 


f 




Quark direct Urea 


d - 


+ u4 


- e _ 


4 P e 






fast 




11 + 


e~ - 


-> d 


4 i/ e 










s — 


-> it 4 




4 v e 










11 + 


e~ - 


— ► s 


4 t^ e 








7r~ condensate 


n+ 


< 7T 


" > 


— > n 


4 e 


~ + «4 


fast 


K~ condensate 


n+ 


< K- > 


— » 71 


4 e +n e 


fast 


Quark bremsstrahlung 


Qi 


+ Q: 


j — > 


Qi4 


Q 2 


+ t/ + t^ 


slow 


Core bremsstrahlung 


n 4 


- n — 


» n - 


-n4 


Ve - 


f V e 


slow 




n4 


-P - 


n + 


-P + 


f e + 


~ V e 






e~ 


+ P - 


-> e~ 


" +P 


4 V 


e + V e 




Crust bremsstrahlung 


e~ 


+ (A,Z) 


— > e~ 


- 4 


{A, Z)+v e - 


- P e slow 



composition. The possible existence of Bose-Einstein condensates or quark mat- 
ter enhance the neutrino emissivity from the core, leading to a more rapid early 
cooling (Table |]) . Superfluidity, on the other hand, has the opposite effect on cool- 
ing. Quantitative constraints have been hampered by the relatively small number 
of young pulsars known, the complication that several of them also display non- 
thermal, beamed X-ray emission from their magnetospheres, and uncertainties in 
distance and interstellar absorption. Table [l^ gives an overview of observed neu- 
tron star temperatures. Because of its unusually low temperature which appears 



to falls well below predictions from standard cooling models£23, the neutron star in 



supernova remnant (SNR) 3C58 received much attention recently (see Sect. I 

4 Properties of Neutron Star Models 

4-1 Nonrotating Star Models 

For many studies of neutron star properties it is sufficient to treat neutron star 
matter as a perfect fluid. The energy-momentum tensor of such a fluid is given by 

=u^ l u v (e + P) + g"" P , (37) 

where m m and u v are four- velocities defined as 

dr^ dr v 
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Table 10. Ages r (in years) and luminosities L (in erg s 1 ) / temperatures T (in K) of pulsars. 
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They are the components of the macroscopic velocity of the stellar matter with 
respect to the actual coordinate system that is being used to derive the stellar 
equilibrium equations. The production of curvature by the star's mass is specified 
by Einstein's field equations, 

G^u = SttT^ , where G^ v = - i g^ v R (39) 

is the Einstein tensor. The Ricci tensor i? Mt , is obtained from the Riemann tensor 
R T iiov by contraction, that is, R^u — R T 9° V which leads to 

with TZ U the Christoffel defined as 



The scalar curvature of spacetime R in Eq. (|39|), also known as Ricci scalar, follows 
from Eq. as 

R = R„ v g^. (42) 

Finally, we need to specify the metric of a non-rotating body in general relativity 
theory. Assuming spherical symmetry, the metric has the form 

ds 2 = „ e 2 *(r) dt 2 + e 2 A(r) ^2 + r 2 ^2 + r 2 ^2 > (43) 

where $(r) and A(r) are radially varying metric functions. Introducing the covari- 
ant components of the metric tensor, 

9tt = ~ e 2 * , g rr = e 2A , g gg = r 2 , g^ = r 2 sin 2 # , (44) 



paper: submitted to World Scientific on February 5, 2008 22 



the non-vanishing Christoffcl symbols of a spherically symmetric body are 

tt c * J x tr Y i i rr 

COS 6* 



s]]lW ^ = -rsin^e^ A , r^ = -siii^ cos^, (45) 

where accents denote differentiation with respect to the radial coordinate. From 
Eqs. (|37l), ( |39l) and (|4^) one derives the structure equations of spherically symmetric 
neutron stars known as Tolman-Oppenheimer-Volkoff equations, EI9uIj 

dP _ e(r)m(r) [1 + P(r)/e(r)] [l + 47rr 3 P(r)/m(r)] 
dr r 2 f — 2m(r)/r 

Note that we use units for which the gravitational constant and velocity of light 
are G = c = 1 so that M Q = 1.5 km. The boundary condition to ([46| ) is P(r = 
0) = P c — P(e c ), where e c denotes the energy density at the star's center, which 
constitutes an input parameter. The pressure is to be computed out to that radial 
distance where P(r = R) = which determines the star's radius R. The mass 
contained in a sphere of radius r (< R), denoted by m(r), follows as 

r 

m(r) =4:7r dr' r' 2 e(r') . (47) 







The star's total gravitational mass is thus given by M = m(R). 

Figure [l4] exhibits the gravitational masses of non-rotating neutron stars as a 
function of central energy density for a sample of equations of state of Tables [l] 
and Each star sequence is shown up to densities that are slightly larger than 
those of the maximum-mass star (indicated by tick marks) of each sequence. Stars 
beyond the mass peak are unstable against radial oscillations and thus cannot exist 
stably (collapse to black holes) in nature. One sees that all equations of state 
are able to support non-rotating neutron star models of gravitational masses M > 
A/(PSR 1913 + 16). On the other hand, rather massive stars of say M ^ 2 M Q can 
only be obtained for those models of the equation of state that exhibit a rather stiff 
behavior at supernuclear densities. The largest maximum mass value, M — 2.2 M@, 
is obtained for HFV, which is caused by the stiffening exchange term. Knowledge 
of the maximum mass value is of great importance for two reasons. Firstly, because 
the largest known neutron star mass imposes a lower bound on the maximum mass 
of a theoretical model. The current lower bound is about 1.60 Mq [neutron star 
4U 0900-40 (= Vela X-l)] which does not set too stringent a constraint on the 
nuclear equation of state. The situation could easily change if an accurate future 
determination of the mass of this neutron star should result in a value that is close 
to its present upper bound of 2.1 M Q . In this case most of the equations of state of 
our collection would be ruled out. The second reaspiijs-that the maximum mass is 



essential in order to identify black hole candidates JHMiS For example, if the mass 
of a compact companion of an optical star is determined to exceed the maximum 
mass of a neutron star it must be a black hole. Since the maximum mass of stable 
neutron stars studied here is ~ 2.2 Mq, compact companions being more massive 
than that value are predicted to be black holes. 
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Figure 14. Non-rotating neutron star mass as 
a function of central density.Q 



Figure 15. Non-rotating neutron star mass as 
a function of redshift. | 



The neutron star mass as a function of gravitational redshift, denned in Eq. fl36|), 
is shown in Fig. [l5|. Maximum-mass stars have redshifts in the range 0.4 < z < 0.8, 
depending on the stiffness of the equation of state. Neutron stars of typically 
M w 1.5 Mq (e.g., PSR 1913+16) are predicted to have redshifts in the considerably 
narrower range 0.2 < z < 0.32. The rectangle covers masses and redshifts in the 
ranges of 1.30 < M/Mq < 1.65 and 0.25 < z < 0.35, respectively. The former range 
has been determined from observational data on X-ray burst source MXB 1636- 
536 L3, while the latter is based on the neutron star redshift data base provided by 
measurements of gamma-ray burst pair annihilation lines .Ej From the redshift value 
of SNR N49 one would expect a neutron mass star of mass 1.1 IS M/Mq IS 1.6. 
Figure [l(] displays the radius as a function of gravitational redshift. 
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Figure [It] shows the moment of inertia of neutron stars, given by 

tt/2 R(8) 


as a function of gravitational mass. In Newtonian mechanics one has for a sphere 
of uniform density / oc R 2 M. The general relativistic expression for the moment 
of inertia, given in Eq. (Eq), is considerably more complicated as it accounts for 



the dragging effect of the local inertial frames (see Fig. 18) and the curvature of 
space-time.cl Nevertheless the qualitative dependence of / on mass and radius as 
expressed in the classical expression remains valid. Estimates for the upper and 
lower bounds on the moment of inertia of the Crab pulsar derived from the pulsar's 
energy loss rate (labeled Rud72), and the lower bound on the moment of inertia 
derived from the luminosity of the Crab nebula (labeled Crab) are shown in Fig. 
0. The arrows refer only to the value of the moment of inertial of the Crab pulsar 
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and not to its mass, which is not known. 



4-2 Rotating Star Models 

The stellar equations describing rotating compact stars are considerably more com- 
plicated than those of non-rotating compact stars These complications have their 
cause in the deformation of rotating stars plus the general relativistic effect pLthc 
dragging of local inertial frames. This reflects itself in a metric of the formBtiJ 

ds 2 = - e 2 »dt 2 + e 2 ^ (dd> -udtf + e 2 »d9 2 + e 2X dr 2 , (49) 

where each metric function v, ip, fx and A depends on the radial coordinate r, 
polar angle 9, and implicitly on the star's angular velocity f2. The quantity u 
denotes the angular velocity of the local inertial frames, which are dragged along 
in the direction of the star's rotation. This frequency too depends on r, 9 and 
f2. Of particular interest is the relative frame dragging frequency Co defined as 
u>(r, 9,0.) = ft — uj(r 7 9, ft), which typically increases from about 15% at the surface 
to about 60% at the center of a neutron star that rotates at its Kepler frequency.H 



4- 3 Mass shedding versus Gravity Wave Emission 

The Kepler frequency, fix, is the maximum frequency a star can have before mass 
loss (mass shedding) at the equator sets in. It sets an absolute upper limit on 
stable rapid rotation. In classical mechanics the expression for fix, determined by 
the equality between centrifuge and gravity, is readily obtained as f^K = \/M/R 3 . 
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Its general relativistic counterpart is given b; 



115 



24? 



V 



h e^~ v 



Pk = 



2n 



(50) 



The primes denote derivatives with respect to the Schwarzschild radial coordinate. 
In order to construct stellar models that rotate at their respective Kepler frequen- 
cies, equation ((5^) is to be solved self-consistently together with Einstein's field 
equations for a given model for the equation of state.oEL 5 



Figure [L9| shows J1k as a function of rotating star mass. The rectangle indicates 
both the approximate range of observed neutron star masses as well as the observed 
rotational periods which, currently, are P > 1.6 ms. One sees that all pulsars so far 
observed rotate below the mass shedding frequency and so can be interpreted as 
rotating neutron stars. Half- millisecond. periods or even smaller ones are excluded 
for neutron stars of mass 1.4 M Q .Hiiia'tia The situation appears to be very different 
for stars made up of self-bound strange quark matter, the so-called strange stars 
which will be introduced in Sect. [|. Such stars can withstand stable rotation against 
mass shedding down to rotational periods in the half-millisecond regime or even 
belowJllJ Consequently, the possible future discovery of a single sub-millisecond 
pulsar spinning at say 0.5 ms could give a strong hint that strange stars actually 
exist, and that the deconfined self-bound phase of 3-flavor strange quark matter is 
in fact the true ground state of the strong interaction rather than nuclear matter. 

The gravitational radiation-reaction driven instability in compact stars sets a 
more stringent limit on rapid stellar rotation than mass shedding, as can be seen 
from Figs. ^0] and [2l]. This instability originates from counter-rotating surface 
vibrational modes, which at sufficiently high rotational star frequencies are dragged 
forward. In this case, gravitational radiation which inevitably accompanies the 
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Figure 19. Onset of mass shedding from rapidly spinning neutron stars, computed for a collection 
of equations of state.p The Kepler period is defined in Eq. (pfj[). 



asphepcai-transport of matter does not damp the instability modes but rather drives 
them.tia'tia Viscosity plays the important role of damping such gravitational- wave 
radiation reaction-driven instabilities at a sufficiently reduced rotational frequency 
such that the viscous damping rate and power in gravity waves are comparable.!!^ 
The instability modes are taken to have the dependence exp[iw m (0)i + irrvf) — 
t/r m (f2)], where uj m is the frequency of the surface mode which depends on the 
angular velocity O of the star, <f> denotes the azimuthal angle, and r m is the time 
scale for the mode which determines its growth or damping. The rotation frequency 
VL at which it changes sign is the critical frequency for the particular mode, m 
(=2,3 p 4_..). It is conveniently expressed as the frequency, denoted by f2^, that 



solves!^ 



or — 



,(0) 



m 



I 3,771 



where 



^(0)../ 2m(m - 1)M 



(51) 



(52) 



2 m + 1 R 3 

is the frequency of the vibrational mode in a non-rotating star. The time scales for 
the gravitational radiation reaction to grow, T Qtmi and for viscous damping to set 
in, T„, m , are given byS 

r a . m = ? ^ m r 1)[ ^ +1 F (J m+ \X R, (53) 



3 (m + l)(m + 2) \2m(m - 1) J \M 
R 2 1 



(2m + 1) (to - 1) v 



(54) 
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Figure 20. Gravitational radiation-reaction 
instability period P T versus mass forrjiewly 
born stars of temperature T = 10 10 K.Q 
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Figure 21. Gravitational radiation-reaction 
instability period P T versus mass for old stars 
of temperature T = 10 6 K.Q 



The shear viscosity, u, depends on a star's temperature according to v{T) oc T -2 . It 
is small in very hot (T « 10 10 K) and therefore young stars but considerably larger 
in cold ones (T sa 10 6 K). Information jabaul the pulsations of the rotating star is 



contained in the functions a m and ^.E^'EEI A characteristic feature of equations 



© through §1) is that merely depends on radius and mass (R and M) of the 
spherical star model, which renders solving these equations rather straightforward. 

Figure ^ shows the critical rotational periods at which the emission of gravity 
waves sets in in hot (T = 10 10 K) pulsars newly born in supernova explosions. Fig- 
ure |H is the analog to Fig. but for old and therefore cold pulsars of temperature 
T = 10 6 K, like neutron stars in binary systems that are being spun up by mass 
accretion from a companion. One sees that the limiting rotational periods P T 
(= 2-k/VI'L) are the smaller the more massive (the smaller the radius) the star 
(cf. Fig. [16]). A comparison between Figs. ^0| and ^l] shows that the instability 
periods are shifted toward smaller values the colder the star due to the larger vis- 
cosity in colder objects. Consequently, the instability modes of neutron stars in 
binary systems are excited at smaller rotational periods than for hot, newly born 
pulsars in supernovae. This conclusion hinges on the role of bulk viscosity of neu- 
tron star matter which may become very large at high temperatures. Sawyer has 
pointed out that it proportional as the sixth poster of the temperature as compared 



with a T~ 2 dependence for the shear viscosityJH3 This means that at temperatures 
T c 10 9 K the bulk viscosity may dominate over the shear viscosity and regulate 
the gravitational radiation-reaction driven instability in rapidly rotating neutron 
stars, pushing their critical rotational periods toward smaller values, possibly even 
as small as the Kepler period. 



The dependence of P T on the equation of state is shown too in Figs. |20| and 21 
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Figure 22. Equation of state of a strange star surrounded by a nuclear crust. Pdrip ( e drip ) denotes 
the pressure at the maximum possible inner crust density determined by neutron drip, e C rust = 
0.24 MeV/fm 3 (4.3 X 10 11 g/cm 3 ). Any inner crust value smaller than that is possible. As an 
example, we show the equation of state for e C rust — 10~ 4 MeV/fm 3 (10 8 g/cm 3 ) " 



too. One sees that the lower limits on P T are set by non-relativistic equations of 
state due to the relatively small radii obtained for stars constructed from them. The 
rectangles in Figs. ^ and |l| labeled "observed" cover both the range of observed 
neutron star masses, 1.1 < M/Mq 5 1.8 as well as observed pulsar periods, P > 1.6 
ms. One sees that even the most rapidly rotating pulsars so far observed have 
rotational periods larger that those at which the gravitational radiation reaction- 
driven instability sets in. Thus all observed pulsars can be understood as rotating 
neutron stars. The observation of pulsars possessing masses in the observed range 
but rotational periods that are smaller than say ~ 1 ms, depending on temperature 
and thus on the pulsar's history, would be in contradiction to our collection of 
equations of state. 

Finally, we mention the r-mode instability which has attracted a great deal 
of attention over the last several years. Like the /-modes discussed above, this 
instability too is produced by the Chandrasekhar, Friedman and Schutz (CFS) 
render every star made of a perfect fluid unstable at all rates 



mechanism 

of rotation l 124 H 125 l The question whether or not it provides a much more severe 
constraint on the stable rotation of viscous stars, however, appears to be still an 
open issue. 

5 Strange Quark Matter Stars 

5.1 The Strange Matter Hypothesis 

The theoretical possibility that strange quark matter may be absolutely stable 
with respect to nuclear matter, that is, the energy per baryon of such matter is 
less tha 
Witten. 



3fl y p y, has been pointed out independently by Bodmcr, Terazawa, and 
This so-called strange matter hypothesis constitutes one of the 
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most startling possibilities regarding the behavior of superdense matter, which, if 
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Figure 23. Relative quark and lepton densiticsr+a cold quark star matter as a function of energy 
density for a bag constant of B 1 / 4 = 145 MeVO 



true, would have implications of fundamental importance for cosmology, the early 
universe, its evolution to the present day, astrophysical compact objects. Even to 
the present day there is no sound scientific basis on which one can either confirm 
or reject the hypothesis so that it-reinains a serious possibility of fundamental 
significance for various phenomena.u'UHj 



5.2 Quark-lepton Composition of Strange Matter 

The relative quark-lepton composition of quark-star matter at zero temperature is 
shown in Fig. |3[ All quark flavor states that become populated at the densities 
shown are taken into account. Strange and charm quark masses of respectively 
0.15 GeV and 1.2 GeV are assumed. Since starsJn their lowest energy state are 
electrically charge neutral to very high precision E£3, any net positive quark charge 
must be balanced by leptons. In general, as can be seen in Fig. there is only 
little need for leptons, since charge neutrality can be achieved essentially among the 
quarks themselves. The concentration of electrons is largest at the lower densities 
of Fig. ^ due to the finite s-quark mass which leads to a deficit of net negative 
quark charge, and at densities beyond which the c-quark state becomes populated 
which increases the net positive quark charge. 

The presence of electrons in strange quark matter is crucial for the possible 
existence of a nuclear crust on such objects, as will be discussed below. Recently 
it has been argued that strange quark matter is a color superconductor (see Sect. 
^) which, at extremely high densities, is in the Color-Flavor-Locked (CF.L1 phase. 



This phase is rigorously electrically neutral with no electrons required.lHIl If the 
CFL phase would extend all the way to the surface of a strange star, then strange 
stars would not be able to carry nuclear crusts because of the missing electric dipolc 
layer. However, for sufficiently large strange quark masses, the "low" density regime 
of strange matter is rather expected to form a 2-flavor color superconductor (2SC) 
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in which electrons are 

As shown in Refs. li3'E£3IH3, the electrons, because they are bound to strange 
matter by the Coulomb force rather than the strong force, extend several hundred 
fermi beyond the surface of the strange star. Associated with this electron displace- 
ment is a electric dipole layer which can support, out of contact with the surface 
of the strange star, a crust of nuclear material, which it polarizes l 134 H 135 l The maxi- 
mal possible density at the base of the crust (inner crust density) is determined by 
neutron drip, which occurs at about 4.3 x 10 11 g/cm 3 . 



5.3 Equation of State of Strange Stars with Crust 



The somewhat complicated situation of the structure of a strange star with crust 
described above can he-represented by a proper choice of equation of state which 



consists of two parts.Eiil At densities below neutron drip it is represented by the 

pal tin cl oar matter, for which we use the 
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The star's strange matter core is 



low-density equation of state of charge-neutr. 
Baym-Pethick-Sutherland equation of state, 
described by the bag model. The graphical illustration of such an equation of state 
is shown in Fig. |2^. Notice that there is a discontinuity in energy density between 
strange quark matter and hadronic matter across the electron surface (dipole gap) 
inside the star where the pressure of the— hadronic crust at its base equals the 



pressure of the strange core at its surfaceJlii 



5.4 Strange and Charm Stars 

Figure ^4] shows the mass of quark matter stars as a function on central star density 
for different, representative bag constants. Stars with central densities smaller than 
a few times 10 15 g/cm 3 are composed of up, down and strange quarks only and are 
referred to as strange stars. Charm quarks are present at central densities higher 
than ~ 10 17 g/cm 3 . Stars with densities larger than that are thus denoted charm 
stars. Only the sequences computed for 5 1 / 4 = 145 MeV consist of matter that 
is absolutely stable at zero external pressure. The other two sequences correspond 
to strange matter whose energy per baryon number is about 960 MeV (i? 1 / 4 = 
167 MeV) and 1030 MeV (B 1/4 = 180 MeV). A stability analysis shows that charm 
stars are—unstable against radial oscillations and thus cannot exist stably in the 
uni verse. Ei^j So the only type of quark stars that could possess physical significance 
are strange stars. 



5. 5 Complete Sequences of Strange Matter Stars 

Since the nuclear crusts surrounding the cores of strange stars are bound by the 
gravitational force rather than confinement, the mass-radius relationship of strange 
matter stars with crusts is qualitatively similar to the one of purely gravitationally 
bound compact stars, that is, neutron stars and white dwarfs. This is illustrated in 
Fig. [2^. The strange star sequence is computed for the maximal possible inner crust 
density, e crustl = edrip- Of course there are other possible sequences of strange stars 
with any smaller value of inner crust density, which will be discussed below. From 
the maximum mass star, the central density decreases monotonically through the 
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Figure 24. Gravitational mass of strange and charm stars for different B 1 / 4 values. 



sequence in each case. The neutron star sequence is computed for HFV of Table [| 
a representative model for the equation of state of neutron star matter, which has 
been combined at subnuclear densities with the Baym-Pethick-Sutherland equation 
of state. Hence the white dwarfs shown in Fig. |2^ are computed for the latter 
equation of state. Those gravitationally bound stars with radii 5 200 km and 
C 3000 km represent stable neutron stars and white dwarfs, respectively. The 
fact that strange stars with crust possess smaller radii than neutron stars leads 
to smaller mass shedding periods, as already indicated by the classical expression 
Pk = 2tt-\/ R 3 /M. Of course the general relativistic expression (|50|), which is to be 

applied to neutron and strange stars, is considerably more complicated However 

the qualitative dependence of Pk on mass and radius remains valid.E£3 Due to 
the smaller radii of strange stars, the complete sequence of such objects (and not 
just those close to the mass peak, as is the case for neutron stars) can sustain 
extremely rapid rotation.tlj In particular, a strange star with a typical pulsar mass 
of ~ 1.45 Mq can rotate at Kepler periods as small as 0.55 5 Px/msec 5 0.8, 

Jl4cj r 



depending on crust thickness and bag constant. I t! " This range-is to be compared 



141 



with Pk ~ 1 msec obtained for neutron stars of the same mass. 

The minimum-mass configuration of the strange star sequence, labeled "a" in 
Fig. |l9|, has a mass of about ~ 0.017 M Q (about 17 Jupiter masses). Depending on 
the chosen value of inner crust densitiz^trange matter stars can be even by orders 



of magnitude lighter than this valuelifj If abundant enough in our Galaxy, such 



low-mass strange stars could be seen by the gravitational microlensing searches. 
Strange stars located to the right of "a" consist of small strange cores (< 3 km) 
surrounded by a thick nuclear crust made up of white dwarf material. Such objects 
are thus called strange dwarfs. Their cores have shrunk to zero at the points labeled 
"X" . What is left is an ordinary white dwarf with a central density equal to the 
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Figure 25. Mass versus radius of strange stars with nuclear crusts (dashed and dotted curves) and 
neutron stars and white dwarfs (solid curve). The strange stars carry nuclear crusts with chosen 
inner densities of e crus t = 4 X 10 11 g/cm 3 and 10 8 g/cm 3 . The two crosses denote the termination 
points of the strange star sequences where the quark matter cores have shrunk to zero. Dote refer 
to maximum mass stars, minimum mass stars are located at the vertical bars labeled "a" £1 



inner crust density of the former strange dwarf .LfjclHI A detailed stability analysis of 
strange stars against radial oscillations shows that the strange dwarfs between "b" 
and "d" are unstable against the fundamental eigenmode.tlj Hence such objects 
cannot exist stably in nature. However all other stars of this sequence are stable 
against oscillations. So, in contrast to neutron stars and white dwarfs, the-brarujies 
of strange stars and strange dwarfs are stably connected with each other. 



14( 



1 13 



Until recently, only rather vague tests of the theoretical mass-radius relation of 
white dwarfs were possible. This has changed dramatically, because of the avail- 
ability of new data emerging from the Hipparcos project.Eda These data allow the 
first accurate measurements of white dwarf distances and, as a result, establishing 
the mass-radius relation of such objects empirically. Figure ^6] shows a comparison 
of several data from the Hipparcos project with the mass-radius relationships of 
strange dwarfs (solid lines) and ordinary white dwarfs computed for different com- 
positions. The outcome could suggests that there is evidence for a bimodality in 
the white-dwarf populations one of which may contain strange-matter cores. 



5.6 Hadronic Crust on Strange Stars and Pulsar Glitches 

Of considerable relevance for the viability of the strange matter hypothesis is the 
question of whether strange stars can exhibit glitches in rotation frequency. From 
the study performed in Refs. 



it is known that the ratio of the crustal moment 
of inertia to the total moment of inertia, / C rust/^totai, varies between 10 -3 and 
~ 10~ 5 . If the angular momentum of the pulsar is conserved in a stellar quake 
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Figure 26. Comparison of theoretical mass-radius relationships of strange dwarfs (e C rust = tdrip) 
and several ordinary white dwarf stars ( 12 C dwarfs:i— dot-dashed, 24 Mg: dotted, 56 Fe: dashed 
double-dotted) with data from the Hipparcos project. £13 



then the relative frequency change and moment oLiujertia change are equal, and 
one arrives for the change of the star's frequency at 
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^ = ^ > ^ ~ (10- 5 -10- 3 )/, with 0</<l. (55) 

Here Iq denotes the moment of inertia of that part of the star whose frequency is 
changed in the quake. It might be that of the crust only, or some fraction, or all of 
the star. The factor / in Eq. (^5|) represents the fraction of the crustal moment of 
inertia that is altered in the quake, i.e., / = |A7|// crust . Since the observed glitches 
have relative frequency changes Ail/57 = (10~ 9 — 1CP 6 ), a change in the crustal 

• 
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moment of inertia of / 5 0.1 would cause a giant glitch even in the least favorable 



case- 



Moreover it turns out that the observed range of the fractional change in 
the spin-down rate, f2, is consistent with the crust having the small moment of 
inertia calculated and the quake invobzing only a small fraction / of that, just as 



in Eq. (pa). For this purpose we writeEli 



where use of Eq. ( p5[) has been made. Equation ( p6|) yields a small / value, i.e., 
/ < (10~ 4 to 10 _1 ), in agreement with / IS 1CP 1 established just above. Here 
measured values of the ratio (Af2/fi)/(AQ/f2) ~ 10~ 6 to 10 -4 for the Crab and 
Vela pulsars, respectively, have been used. 
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Figure 27. Cooling behavior of a 1.4 Mq neutron star based on competing assumptions about 
the behavior of superdense matter. Three distinct cooling scena rios , referred to as "standard", 
"intermediate", and "enhanced" can be distinguished (cf. Sect. p.7[ ). The barrnl-like structures 
reflects the uncertainties inherent in the equation of state of superdense matter 



6 Thermal evolution of Neutron and Strange Stars 

The predominant cooling mechanism of hot (temperatures of several ~ 10 10 K) 
newly formed neutron stars immediately after formation is neutrino emission, with 
an initial cooling time scale of seconds. Already a few minutes after birth, the 
internal neutron star temperature drops to ~ 10 9 K. Photon emission overtakes 
neutrinos only when the internal temperature has fallen to ~ 10 8 K, with a corre- 
sponding surface temperature roughly two orders of magnitude smaller. Neutrino 
cooling dominates for at least the first 10 3 years, and typically for much longer in 
standard cooling (modified Urea) calculations. Being sensitive to the adopted nu- 
clear equation of state, the neutron star mass, the assumed magnetic field strength, 
the possible existence of superfluidity, meson condensates and quark matter, theo- 
retical cooling calculations, as summarized in Fig. |27], provide most valuable infor- 
mation about the interior matter and neutron star structure. The thermal evolution 
of a neutron star also yields information about such temperature sensitive prop- 
erties as transport coefficients, transition to superfluid states, crust solidification, 
and internal pulsar healing mechanisms such as frictional dissipation at the crust- 
superfluid interfaces .E13 The stellar cooling tracks in Fig. p7] are computed for the 
broad collection of equations of state introduced in Sect.gy Recent Chandra ob- 
servations have identified pulsar J0205+6449 at the center of the young Crab-like 
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supernova remnant 3C58. Historical evidence suggests an association of the rem- 
nant with supernova SN 1181, which makes 3C58 younger (see Table 0) than Crab. 
The temperature of J0205+6449 was recently determined to be log 10 T — 6.0334 
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en from Fig. 27, falls well below predictions of standard cool- 
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15( ], to the direct Urea 



(see Table | for an overview of neutron star cooling mechanisms). 



7 Color Superconductivity of Quark Matter 

There has been much recent progress in our understanding of quark matter, culmi- 



nating in.the Hi.sroy.ery that if quark matter exists it will be in a color supe 



ing state.t22 
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The phase diagram of such matter is very complex 
asymptotic densities the ground state of QCD with a vanishing strange quark mass 
is the color-flavor locked (CFL) phase. This phase is electrically neutraLip bulk 



for a significant range of chemical potentials and strange quark masses.ELIJ If the 



strange quark mass is heavy enough to be ignored, then up and down quarks may 
pair in the two-flavor superconducting (2SC) phase. Other possible condensation 



patters are the CFL-AT phase E£a and the color-spin locked (2SC+s) phaseMS The 



magnitude of the gap energy lies between ~ 50 and 100 MeV. Color superconduc- 
tivity thus modifies the equation of state at the order (A//i) 2 level, which is only a 
few percent. Such small effects can be safely neglected in present determinations of 
models for the equation of state of neutron stars and strange quark matter stars, as 
is the case here. There has been much recent work on howxaJpr h" 1 HP rf, L nrl 
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activity 
15 i and 



in neutron stars could affect their properties. (See Refs 
references therein.) These studies reveal that possible signatures include the cool- 
ing by neutrino emission, the pattern of the arrival times of supernova neutrinos, 
the evolution of neutron star magnetic fields, rotational (r-mode) instabilities, and 
glitches in rotation frequencies. 

Aside from neutron star properties, an additional test of color superconductivity 
may hep-provided by upcoming-cosmic ray space experiments such as AMS E23 and 
ECCO.ESJ As shown in Ref. Ii£3, finite lumps of color-flavor locked strange quark 
matter (strangelets) , which should be present in cosmic rays if strange matter is the 
ground state of the strong interaction, turn out to be significantly more stable than 
strangelets without color-flavor locking for wide ranges of parameters. In addition, 
strangelets made of CFL strange matter obey a charge- mass relation of Z/A oc 
A -1 / 3 , which differs significantly from the charge-mass relation of strangelets made 

ild, be constant for 
This .difference 



of "ordinary" strange quark matter. In the latter case, ZlA 
small baryon numbers A and Z/A cx A~ 2 / 3 for large AL 
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L6S 



may allow an experimental test of CFL locking in strange quark matter Jl£3 



8 Particle Thresholds and Quark Deconfinement in Rotating 
Neutron Stars 

Figures ^ and |2^ reveal that the weakening of centrifuge accompanied by the 
slowing-down of a rotating neutron star causes a significant increase of its central 
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Figure 28. Central density versus rotational fre- 
quency of neutron stars of constant baryon num- 
ber, A, computed for Gg^OO Threshold den- 
sities of various possible phases of superdense 
matter are indicated. M(0) is the non-rotati»« 
star mass, Op; stands for the Kepler frequency.! 32 ! 



Figure 29. Central star density versus rotational 
frequency for two sample stars computed for 
equation of state Gg 



density, which is accompanied by a rearrangement of the particle population and 
the creation of new phases of matter. From Fig. 29, for instance, one reads off that 
the central density of a neutron star of mass M — 1.42 M@ increases from about 
450 MeV/fm 3 for rotation at fix to more than 1500 MeV/fm 3 for = 0, which is a 
~ 66% effect. Of course, this effect is smaller for the lighter stars, since for them the 
changes of the gravitational pull are weaker. Such dramatic changes in the interior 
density of slowing-down neutron stars imply profound changes of their interior 
composition, because of the changing particle composition with density. Figures 
|30| through 33 illustrate how these changes carry over to the internal structure of 
conventional neutron stars. In each case the star's rotational frequency covers the 
maximal possible range < < Ok- It is evident that as rotating stars spin down 
they become significantly less deformed, which leads to closer equality between the 
polar and equatorial radii. At the same time the central density rises from below to 
above the threshold densities of heavier baryons (i.e., S, A, 5). For some pulsars the 
mass and initial rotational frequency may be such that the central density rises from 
below the critical density for dissolution of baryons into their quark constituents. 
Examples of which are shown in Figs. |34| through [37|. The evolution of the central 
mass densities of these stellar models can be inferred in reference with Fig. ^9|. For 
rotational frequencies below ~ 1250 s" 1 the quark- hybrid star models consist of an 
inner region of pure quark matter (q) surrounded by a few kilometer thick shell 
of mixed phase of hadronic and quark matter. This phase consists of structures 
like hadronic drops (hd), hadronic rods (hr), hadronic slabs (hs), quark slabs (qaL 
quark rods (qr), and quark drops (qd), in each case arranged in a lattice structure. Eil 
This shell is surrounded by a sphere of hadronic liquid (hi), itself enveloped in a 
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thin nuclear crust of heavy ions. These quark-hadron structures are a consequence 
of the competition of the Coulomb and surface energies of the hadronic and quark 
matter phaseH They may have dramatic effects on pulsar observables including 
transport properties and the theory of glitches. 

If the mass and initial rotational frequency of a pulsar is such that during 
its slowing-down phase the interior density rises from below to above the critical 
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density for the quark-hadron phase transition, first at the center where the density 
is highest (Figs. |28| and |2^) and then in a region expanding in the radial outward 
direction away from the star's center, matter will be converted from the relatively 
incompressible nuclear matter phase to the highly compressible quark matter phase. 
The tremendous weight of the overlaying layers of nuclear matter tends to compress 
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the quark matter core, which causes the entire star to shrink on a macroscopic length 
scale, as shown in Figs. |34| through [37| The mass concentration in the core will 
be further enhanced by the increasing gravitational attraction of the quark core on 
the overlaying nuclear matter. The moment of inertia thus decreases anomalously 
with decreasing rotational frequency as the new phase slowly engulfs a growing 
fraction of the star li£j, as can be seem from Figs. ^ and [39[ F igure [58] shows the 
moment of inertia, /, computed self-consistently from Eq. (f48|) for severaLsample 
stars having the same baryon number but different internal constitutions .11^3 The 
curve labeled M — 1.420 M©, computed for Gg^gQ, shows the moment of inertia 
of the quark-hybrid star of Figs. |34| and |35| . The other curves correspond to a 
standard hyperon star constructed for Gjjj^f and a standard neutron star where 
hyperons were ignored. In accordance with what has been said just above, the 
shrinkage of quark-hybrid stars driven by the development of quark matter cores is 
the less pronounced the smaller the quark matter cores which are being built up in 
their centers during spin-down. Correspondingly the reduction of / weakens with 
decreasing star mass, as shown in Fig. ^8| for several sample star masses. 

The decrease of the moment of inertia caused by the quark-hadron phase tran- 
sition shown in Fig. |3^ is superimposed on the response of the stellar shape to a 
decreasing centrifugal force as the star spins down due to the loss of rotational 
energy. In order to conserve angular momentum not carried off by particle radi- 
ation from the star, the deceleration rate ft (< 0) must respond correspondingly 
by decreasing in absolute magnitude. More than that, Q may even change sign, as 



shown in Figs. |39 I 164 ! which carries the important astrophysical information that 
an isolated pulsar may spin up during a certain period of its stellar evolution. Such 
an anomalous decrease of I is analogous to the backbending phenomenon known 
from nuclear physics, in which case the moment of inertia of an atomic nucleus 
changes anomalously because of a change in phase from a nucleon spin aligned 
state at high ang ular mpm |Cntum to a pair correlated supcrfluid phase at low an- 



gular momentum.^ —'i^i The stellar backbending shown in Fig. ^ shows that 
stars evolving from b — > a are rotationally accelerated (f2 > 0), while stars evolving 
from a —> b are rotationally decelerated (fi < 0). The structure in the moment 
of inertia and, specifically, the backbending phenomenon may dramatically modify 
the timing structure of pulsar spin-down, rendering the observation of quark matter 
in neutron stars accessible to radio astronomy. 



9 Evolution of a Pulsar's Braking Index 

Pulsars are identified by their periodic signal believed to be due to a strong magnetic 
field fixed in the star and oriented at an angle from the rotation axis. The star's 
angular velocity of rotation decreases slowly but measurably over time, and usually 
the first and occasionally second time derivative can also be measured. Various 
energy loss mechanisms could be at play such as the magnetic dipole radiation, 
part of which is detected on each revolution, as well as other losses such as ejection 



of charged particles .E£3 Without loss of generality, one can assume that pulsar slow- 



down is governed by a single mechanism, or several mechanisms having the same 
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power law. The energy balance equation can then be written in the form 
dE d (\ 



,// ,// K2 W)*) = -ctr+ (57) 

In the case of magnetic dipole radiation, the constant C is equal to C = |/z 2 sin 2 a 
(fj, denotes the star's magnetic moment), and n = 3 if I is kept constant during 
spin-up (down). If, as is customary, the star's angular velocity f2 is regarded as the 
only time dependent quantity, one obtains the usual formula for the rate of change 
of the pulsar frequency, given by 

Q = — K Q, n , (58) 

with K a constant and n the braking index. As known from Figs. |38| and the 
moment of inertia is not constant in time but responds to changes in rotational 
frequency. This response changes the value of the braking index in a frequency 
dependent manner, that is n — n(Q). Thus during any epoch of observation, the 
braking index will be measured to be different from it canonical value n = 3 by 
a certain amount. How different depends, for any given pulsar, on its rotational 
frequency and for different pulsars of the same frequency, on their mass and on 
their internal constitution. When the frequency response of the moment of inertia 
is taken into account, Eq. (p8|) is to be replaced witbEK 



On I dl_ I dl 

" 1 - do. n i I da 



n = -2iK nr — = - K n n [ i-ffin+{ fin] -•••), (59) 



2I+jL n \ 21 \2I 



where K = C/I. This explicitly shows that the frequency dependence of f2 corre- 
sponding to any mechanism that absorbs (or deposits) rotational energy cannot be 
a simple power law as given in equation (p8[) (with K a constant). It must depend 
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on the mass and internal constitution of the star through the response of the mo- 
ment of inertia to rotation as expressed in ([59|). It is interesting to note that Eq. 
( j59|) can be represented in the form of ( ^8|) but now with a frequency dependent 
prefactor, by evaluating 



~~dT 

dn 



n(Q) = = 3 — ZM f " . (60) 

fJ 2 21+ 



One sees that this effective braking index depends explicitly and implicitly on Q 
and reduces to the canonical expression n = 3 only if / is independent of frequency. 

As an example, we show in Fig. ^ the change of the braking index with fre- 
quency for two selected quark-hybrid stars. For illustration purposes we assume 
dipole radiation. Because of the structure in the moment of inertia, driven by the 
phase transition of baryonic matter into deconfined quark matter, the braking index 
deviates dramatically from n = 3 at those rotational frequencies where quark de- 
confinement leads to the built-up of pure quark matter cores in the centers of these 
stars. Such anomalies in n(n) are not obtained for conventional neutron stars or 
hyperon stars because their moments of inertia increase smoothly with f2, as known 
from Fig. |38|. The observation of such an anomaly in the timing structure of pulsar 
spin-down could thus be interpreted as a possible signal of quark deconfinement in 
the center of a pulsar. Of course, because of the extremely small temporal change 
of a pulsar's rotational period, one cannot measure the shape of the curve which is 
in fact not necessary. Just a single anomalous value of n that differed significantly 
from the canonical value of n — 3 would suffice l+HI 170 

As a final but very important point on this subject, we estimate the typical 
duration over which the braking index is anomalous if deconfinement is well pro- 
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Figure 41. Braking index versus time for the quark- hybrid star of mass M 
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nounced. The time span can be estimated from 

Aft _ AP 
ft ~ P 



AT 



(61) 



where Aft is the frequency interval of the anomaly. The range over which n(ft) 
is smaller than zero and larger than six (Fig. [lo]) is Aft « — 100 s" 1 , or AP 
— 27rAft/ft 2 » 3 x 10~ 4 s at ft = 1370 s _1 . So, for a millisecond pulsar whose 
period derivative is typically P ~ 10 -19 one has AT ~ 10 8 years, as illustrated in 
Fig. fll]. The dipole age of such pulsars is about 10 9 years. So, as a rough estimate 
one may expect about 10% of the 25 presently known solitary millisecond pulsars 
to be in the transition epoch during which pure quark matter cores are gradually 
being built up in their centers. These pulsars could be signaling the ongoing process 
of quark deconfinement in their cores. Last but not least we note that the spin- up 
time (region b-a in Fig. |3^) is about 1/5 of the time span AT, or about 1/50 of the 
dipole age t = — ft/ (ft(n — 1)). 

It is probably needless to say that the observation of a pulsar with an anoma- 
lously large braking index, signaling the existence of quark matter in the centers of 
pulsars, would be a momentous discovery with most far-reaching consequences for 
nuclear and particle physics. It would help to clarify how quark matter behaves, 
and give a boost to theories about the early Universe as well as laboratory searches 
for the production of quark matter in heavy-ion collisions, which has become a 
forefront area of modern physics. On an even more fundamental level, such a dis- 
covery would prove that the essentially free quark state predicted for matter at 
very high energy densities actually exists, and give us a picture of an early phase of 
the Universe that is based on radio pulsar observation, which may be coined quark 
astronomy. 
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Figure 42. Evolution of spin frequencies of acciiettrig X-ray neutron stars with (solid curves) and 
without (dashed curves) quark deconfinement .EHI The spin plateau around 200 Hz signals the 
ongoing process of quark re-confinement in the stellar centers 



10 Quark Matter in Neutron Stars in Low-Mass X-ray Binaries 

The signal of quark deconfinement described in Sect. || applies only to isolated 
neutron stars, where deconfinement is driven by the gradual stellar contraction as 
the star spins down. The situation is reversed in neutron stars in binary systems, 
which experience a spin-up torque because of the transfer of angular momentum 
carried by the matter picked up by the star's magnetic field from the surrounding 
accretion disk. The spin-Hp, torque causes a change in the stars' angular momentum 
according to the relation 



171 



^ = Ml(r m ) - N(r c ) , (62) 



where M denotes the accretion rate and 



!(r m ) = VM^n (63) 

is the angular momentum added to the star per unit mass of accreted matter. The 
quantity N stands for the magnetic plus viscous torque term, 

N(r c ) = K [i 2 r- 3 , (64) 

with /i = B?B the star's magnetic moment. The quantities r m and r c denote the 
radius of the inner edge of the accretion disk and the co- rotating radius, respectively, 
and are given by (£ ~ 1) 

r m = £r A , (65) 
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Figure 43. Calculated spin distribution of X-ray neutron stars. The spike in the calculated distri- 
bution (unshaded diagram) corresponds to the spinout of the quark mattcpflhase. Otherwise the 
spike would be absent. The shaded histogram displays the observed data. £13 



and 

r c = (Mil~ 2 ) 1/3 . (66) 

Accretion will be inhibited by a centrifugal barrier if the neutron star's magne- 
tosphere rotates faster than the Kepler frequency at the magnetosphere. Hence 
T m < T c , otherwise accretion onto the star will cease. The Alfen radius ta in Eq. 
( |65| ) is defined by 

- = (i4i)" 7 - (67 » 

The rate of change of a star's angular frequency f2 (= J /I) then follows from Eq. 
© as 

with the explicit time dependences as indicated. There are two terms on the right- 
hand-side of Eq. (^) that grow linearly respectively quadratically with fi. Ignoring 
the linear term shows that mass transfer can spin up a neutron star to an equilibrium 



period (P = 2tt/^) o£Z3 



-3/7 

M \ ( M 



-5/7 



P "' = 2 - 4 ™ B ^J & ^' ,Br ' <m) 

where Rq and Bg are the star's radius and its magnetic field in units of 10 6 cm 
and 10 9 G, respectively. MEdd m Eq. ( |69[ ) denotes the maximum possible accretion 
rate, known as the Eddington limit, at which the accretion luminosity equals the 
luminosity at which the radiation pressure force on ionized hydrogen plasma near 
the star balances the gravitational acceleration force exerted by the star on the 
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plasma. This condition leads to an Eddington accretion rate of MEdd — 1-5 x 
10 _8 i? 6 Af Q /yr. For a typical accretion rate of M_ 10 = Af /(10~ 10 M o /yr), the 
Eddington rate can be expressed as MEdd = 150 Rq Mz\q M. The low-mass X-ray 
binaries (LMXBs) observed with the Rossi X-ray Timing Explorer are divided into 
Z sources and A(toll) sources, and appear to accrete at rates of M-%o ~ 200 and 
M_io ~ 2, respectively. 

The temporal change of the moment of inertia of accreting neutron stars which 
undergo phase transitions is crucial for their spin evolution.EIil This temporal 
change, on the other hand, renders the computation of the moment of inertia, 
defined in Eq. (ff8|), very cumbersome since each quantity on the right-hand-side 
of ( [48| ) varies accordingly during stellar spin-up. The solution of Eq. ( |68| ) in com- 
bination with Eq. ( p!8| ) is shown in Fig. E^. The result is most striking. One sees 
that quark matter remains relatively dormant in the stellar core until the star has 
been spun up to frequencies at which the central density is about to drop below 
the threshold density at which quark matter exists. As known from Fig. [59[ this 
manifests itself in a significant increase of the star's moment of inertia. The angular 
momentum added to a neutron star during this phase of evolution is therefore con- 
sumed by the star's expansion, inhibiting a further spin-up until the quark matter 
has been converted into a mixed phase of matter made up of hadrons and quarks. 
Such accreters, therefore, tend to spend a greater length of time in the critical 
frequencies than otherwise. There will be an anomalous number of accreters that 
appear at or near the same frequency, as shown in Fig. |4^. This is what was ac- 
tually found recently with the Rossi X-ray Timing Explorer (shown by the shaded 
area in Fig. B^ n Xkiaxk deconfinement constitutes a most striking explanation, ^ n 



this anomaly l 1Y 4[ L73 H 174 l though alternative explanations were suggested too. 
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